Physical Chemistry Chemical Physics

The repair of photo-induced DNA lesions through nucleotide excision repair machinery is still the source of important questions. It has been observed that the repair rate of the different cyclobutane pyrimidine dimers, i.e. the photoproducts induced by dimerization of two {pi}-stacked pyrimidines (T<>T, T<>C, C<>T, C<>C), depends on the nucleobases involved in the lesion. TT derivatives (T<>T) are removed more slowly than those containing cytosine, especially in 5. Using all-atom molecular dynamics simulations and free-energy calculations, we demonstrate that the variation of the repair rate observed in human skin and in cultured cutaneous cell is associated to the recognition of the four lesions by the DDB2 protein moiety, and more specifically by the differential structural deformation induced on the complementary strand. Indeed, while C<>C and C<>T induce a larger deviation on the groove parameters, T<>T and T<>C, instead, affect DNA structure to a lesser extent. less affected. These effects then hamper differentially the downstream recruitment of the repair complexes. The observed DNA deformation correlates with the experimental repair rate and provides a structural rationale for the different repair rates of CPD by nucleotide excision repair machinery. GRAPHICAL ABSTRACT O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=105 SRC="FIGDIR/small/724087v1_ufig1.gif" ALT="Figure 1"> View larger version (43K): org.highwire.dtl.DTLVardef@cf6b6dorg.highwire.dtl.DTLVardef@195e35forg.highwire.dtl.DTLVardef@1829296org.highwire.dtl.DTLVardef@165baba_HPS_FORMAT_FIGEXP M_FIG C_FIG

Photoactivated intermolecular electron transfer (ET) in cryptochromes proceeds along chains of aromatic residues and creates a spatially separated pair of radical electrons. Ultrafast time-dependent spectroscopy can provide experimental insight into this process and theoretical estimates of charge transfer rates are commonly obtained via Marcus theory. Here, we present a new perspective on the ET in European robin cryptochrome 4a (ErCry4a) that synthesizes insights from real-time ET calculations, ultrafast spectroscopic measurements and analytical derivations. The simulations exemplify that molecular vibrations play an essential role in enabling the ET dynamics, which was further rationalized through analytical derivations. Ultrafast pump-probe spectroscopy provided experimental access to the first 1.5 ns of the ET cascade, where multiple radical pair recombination rates arise due to the dynamic equilibrium along the ET chain. We show that the motions of the protein environment and the ET dynamics are inseparably coupled, violating the timescale separation required for Marcus theory. The presented results highlight that non-equilibrium coupling between electronic and nuclear motion dominates ET kinetics in ErCry4a during the first nanosecond after photo-excitation. The findings exemplify the limits of Marcus theory and refine the interpretation of ultrafast spectroscopic signatures in cryptochromes.

In mammalian cells, lipid monolayers support the integrity of lipid droplets (LDs), organelles that function as storage for neutral lipids. Liver-targeting illnesses such as liver cancer interrupt normal LD metabolism and prompt changes in the chemical content of these organelles, which can have effects on structural and organizational behavior of the lipids. In LDs, liver cancer induces concentric crystalline phases of cholesteryl esters (CEs) and triglycerides near the NL-monolayer interface, which become more pronounced as CE concentration increases. Yet, there is little known about how this phenomenon may link to persistence of undigested LDs in liver cancer patients. To shed light on this, all-atom molecular dynamics simulations were used to model LD micropipette aspiration experiments and gain insight into the effect of CE concentration on partitioning, structural, and mechanical properties of LDs. We successfully model micropipette aspiration by application of constant surface tension laterally, which stretches lipid bilayers and monolayers as the magnitude increased. The results show increased phospholipid packing due to insertion of CE fatty tails into the monolayer. Increasing CE concentration induces a non-linear change in surface packing defects on the LDs, notable rigidification, and stiffness. Taken together, these insights improve our understanding of the physical properties at the LD monolayer-core interface during liver cancer progression.

Stereochemistry underlies structure-function relationships across biology and materials science, ranging from proteins to electronic and spintronic materials. In this work, we investigate the electron transport properties of different oligopeptide stereoisomers using experiments and computational modeling. Single-molecule electronic experiments show that stereochemical modifications in tyrosine-based peptides lead to significant variations in molecular conductance along the peptide backbone due to enhanced stacking interactions and electronic coupling of aromatic side chains. In addition, stereochemical variations in alanine-based peptides give rise to changes in conductivity due to secondary structure interactions arising from {beta}-turn conformations. All-atom molecular dynamics (MD) simulations and quantum mechanical calculations are used to understand the molecular origins of the effect of stereochemistry on the structural and electronic properties of peptides. Overall, this work shows that stereochemical modification of non-terminal amino acids effectively controls electron transport due to aromatic side chain interactions or secondary structure effects. These insights open new avenues for the molecular design of peptide-based electronic materials with enhanced function.

Plant-derived extracellular vesicles (PDEVs) have emerged as superior candidates for oral drug delivery, exhibiting a gastrointestinal survivability that significantly exceeds that of mammalian exosomes or synthetic liposomes. However, the biophysical rules governing how plant genomic regulation translates into this exceptional mechanical resilience remain unknown. Here, we present a predictive multiscale model of plant-derived extracellular vesicles, linking a parameterised genetic state space to emergent mesoscale mechanics via supra-molecular coarse-grained molecular dynamics (SCG-MD). We demonstrate that the upregulation of sterol methyltransferases (SMT) during the plants theoretical Defence state drives the formation of a phase-separated composite architecture, where rigid domains occupying approximately 36% of the membrane surface area effectively arrest crack propagation. This state achieves a critical rupture tension of 367.0 {+/-} 0.7 mN m-1 corresponding to a 39% increase over the wild-type Ripening state. Crucially, we find that chemical composition alone is insufficient for this reinforcement; vesicles with actively sorted lipid domains (Seeded topology) outperform randomised mixtures (Spontaneous topology) by 23% at identical concentrations. Furthermore, while fluid vesicles stiffen reactively under gastric acid shock (pH 2.5) due to the steric jamming of thermodynamically neutralised headgroups, the Defence state exhibits mechanical homeostasis. These findings suggest that PDEVs function as genetically tunable composite materials, offering a design blueprint for next-generation bio-inspired drug delivery vectors. Ultimately, these theoretical indices provide a predictive biophysical framework awaiting empirical confirmation via in vitro nanomechanical assays.

The molten globule (MG) state is an intermediate in the unfolding pathway of proteins, typically triggered by denaturing agents such as urea, extreme pH, high pressure, or heat. The microscopic details of such states are far from understood. Here we study the MG states in protein Hen Egg-White Lysozyme (PDB ID: 1AKI) using microscopic constant pH molecular dynamics (CpHMD) simulations and experiments across a wide pH range. We observe that the titratable residues act as key drivers of conformational fluctuations, promoting the emergence of MG states at extreme pH. These states display partial unfolding, and small global structural changes (< 7% deviation). Hydration around the fluctuating acidic residues shows reduced water density and weakened hydrogen bonding at low pH. At high pH, hydration around acidic residues increases relative to pH = 7, whereas hydration around basic residues decreases. The translational and rotational dynamics of hydration water also exhibit pronounced pH dependence: the translational diffusion coefficient (Dtrans) increases linearly with decrease in pH in acidic medium and increases linearly with increasing pH in the basic regime. The rotational diffusion (Drot) shows similar dependencies on pH except a break at pH {approx} 4 corresponding to acidic residue pKa values. Our results may be useful to identify ligand binding of lysozyme in extreme pH conditions.

Spiral geometries commonly occur in natural and engineered systems and are fundamentally described by curvature and torsion. In deformation-dominated systems, these variables evolve dynamically, requiring a continuum mechanical framework to link geometry and deformation. This study focused on refractile bodies (R-bodies), protein supramolecular assemblies that undergo reversible roll-spiral transformations in response to stimuli such as pH changes. Although multiple R-body types with distinct morphologies and unrolling behaviours were experimentally identified, their deformation mechanisms lack quantitative theoretical descriptions. We proposed a deformation-gradient-tensor-based continuum model incorporating geometrical mapping from the rolled to spiral state within a unified framework. The model successfully reconstructed macroscopic deformation behaviours of types 51, 7, and Pa R-bodies by capturing differences in unrolling behaviours, tapered geometry, and spatio-temporal evolution. The analysis revealed that deformation proceeds through a coupled process in which the curvature decreases via straightening, while the torsion increases by twisting. Importantly, the framework connected the macroscopic morphology with microscopic lattice deformation, enabling quantitative inference of lattice intervals and angles. The proposed comprehensive geometric model of the R-body roll-spiral transformation offers a general mathematical foundation for understanding deformation-driven spiral transformations in soft matter systems.

DNA topology is a key regulator of chromatin structure and transcription, yet its direct role in transcription factor recognition remains unclear. Here, we investigate how distinct DNA topological states modulate binding of the Saccharomyces cerevisiae bZIP transcription factor GCN4 using topologically defined plasmids. By combining, complementary biochemical approaches, including Bio-Layer Interferometry applied here for the first time to topology-dependent protein-DNA interactions, we show that DNA supercoiling directly reshapes GCN4-DNA recognition. Positively supercoiled DNA forms more stable and persistent complexes, whereas negatively supercoiled DNA retains greater conformational heterogeneity. To interpret these effects, we performed multiscale molecular simulations. Coarse-grained simulations of plasmids recapitulate the global topology-dependent trends observed experimentally, while matched minicircle models reproduce the same behaviour at the local scale. In strong agreement with experimental data, simulations reveal that DNA topology modulates the conformational ensemble of the GCN4 basic region. Overall, positively supercoiled DNA promotes a more ordered binding mode and localized protein distribution, whereas negatively supercoiled DNA supports increased structural plasticity. These findings identify DNA topology as an active determinant of transcription factor recognition and provide a multiscale framework linking global DNA mechanics to local protein-DNA interactions. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=113 SRC="FIGDIR/small/722604v1_ufig1.gif" ALT="Figure 1"> View larger version (51K): org.highwire.dtl.DTLVardef@18f8ba9org.highwire.dtl.DTLVardef@11a395dorg.highwire.dtl.DTLVardef@ac093borg.highwire.dtl.DTLVardef@923212_HPS_FORMAT_FIGEXP M_FIG C_FIG

HrpZ2, a harpin protein produced by Pseudomonas syringae, a gram-negative plant pathogenic bacterium, elicits hypersensitive response and pathogen defense in non-host plants. Harpins from various bacterial sources elicit varying responses in different non-host plants, due to its structural variations, their precise mechanisms of action are not yet completely understood. As per previous reports, harpins from diverse bacterial sources interact with distinctive members of integral membrane proteins, known as aquaporins. For example, harpin (Hpa1Xoo) interacts with OsPIP1;3 in rice, whereas, in Arabidopsis the harpins Hpa1 and HrpZ interacts with AtPIP1;4 and AtPIP1;3 respectively. Here, we conducted the first genome-wide computational screening of protein-protein interactions between HrpZ2 and all 47 members of tomato aquaporins. Molecular docking identified nine interactors across five subfamilies of aquaporins, with HrpZ2 N-terminal residues mediating these interactions. We validated these via molecular dynamics (MD) simulations, principal component analysis, and free energy landscape analysis, assessing the stability (RMSD, RMSF, radius of gyration), dynamics, and affinity (MM-GMSA). PIP complexes, especially PIP2;1 (-460.46 kcal/mol) and PIP1;7 (-303.82 kcal/mol), exhibited superior stability, compactness, and defined energy minima, confirming PIPs as primary sensors of harpins. Non-PIP aquaporins like TIP1;1 and NIP4;1 showed moderate stability, outperforming weaker interactors (SIP2;1, XIP1;5, XIP1;3). These findings provide robust evidence that HrpZ2 preferentially targets PIPs in tomato, while engaging TIPs and NIPs as auxiliary partners. This multifaceted interaction profile of harpins suggests complex plant-pathogen recognition, modulating aquaporin-mediated cellular responses like growth and stress management in plants.

HIV-1 buds from infected cells as immature virion particles with a scattered envelope glycoprotein (Env) distribution on their envelope. It then undergoes maturation, during which the viral protease cleaves the Gag polyprotein at multiple sites, leading to structural reorganization of the viral particle and lateral redistribution of Env proteins, ultimately rendering the virion infectious. However, the underlying mechanism of maturation-induced Env reorganization remains elusive. In this study, we combine microsecond-long all-atom (AA), bottom-up coarse-grained (CG) molecular dynamics simulations, and diffusion model-based backmapping to investigate the structural organization and key interactions of Env in viral membranes. AA simulations of fully glycosylated Env embedded in HIV-1 mimetic asymmetric bilayers were first performed to characterize its conformational dynamics and Env-lipid interactions. We then developed a bottom-up CG model of glycosylated Env from that AA data and simulated the mature HIV-1 virion envelope containing multiple Env proteins. The CG simulations predict that Env proteins form clusters through interactions mediated by the cytoplasmic tail domain (CTD) and adopt diverse tilted conformations within these clusters. These CG simulations were then backmapped to AA resolution and further AA simulations were carried out to identify, in detail, the specific interacting residues in the Env clusters. Additionally, analysis of epitope accessibility shows that broadly neutralizing antibodies (bnAbs) targeting the V1/V2 and V3 loops may efficiently interact with Env clusters on the mature virion surface. Together, these results provide a molecular mechanism for Env oligomerization during viral maturation and offer new insights into the accessibility of bnAb epitopes on Env clusters.

Biomolecular condensates (BMCs) formed through liquid-liquid phase separation (LLPS) play key roles in cellular organization, yet their molecular-scale interactions with membranes remain poorly understood. Here, we use surface-enhanced infrared absorption spectroscopy (SEIRAS) to quantify the secondary structure of poly-GR (glycine-arginine) inside condensates in contact supported lipid bilayers. Amide I SEIRAS provides a surface-sensitive measure of the peptide backbone conformation at the membrane-condensate interface. Results show that neutral POPC bilayers preserve the {beta}-sheet rich organization of the bulk condensate, whereas negatively charged POPC/POPS bilayers reduce {beta}-sheet content, and enhanced {beta}-turn formation, accompanied by perturbations of the lipid headgroups. These results demonstrate that membrane charge modulates condensate secondary structure, hydration, and interfacial behavior, providing molecular-level insight into electrostatic regulation of condensate organization at the membrane-water interface. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=199 SRC="FIGDIR/small/726873v1_ufig1.gif" ALT="Figure 1"> View larger version (112K): org.highwire.dtl.DTLVardef@9d9e27org.highwire.dtl.DTLVardef@1fa3fe4org.highwire.dtl.DTLVardef@78bcceorg.highwire.dtl.DTLVardef@b3a2ae_HPS_FORMAT_FIGEXP M_FIG C_FIG

The isolation of cellulose nanofibrils (CNFs), a promising precursor for sustainable and high-performance materials, has relied on chemically intensive, energy-demanding processes. As these processes were originally designed for the isolation of CNFs from wood, we herein show that the intrinsic ultrastructure of non-structural plant cells provides unique opportunities, namely direct access to loosely organized cellulose nanonetworks. We demonstrate that this loose nanofibrillar tissue can be transformed into CNFs with sizes down to elementary nanofibrils ([~]4 nm) at high yields (reaching [~]32%) under exceptionally mild hydrothermal conditions. Three distinct plants were evaluated and the physicochemical properties of the obtained nanonetworks and corresponding CNFs were thoroughly studied, including the hydrodynamics of the resulting gels. Films prepared from the obtained CNFs showed similar performance to those obtained from conventionally isolated wood-based CNFs. Overall, this study demonstrates that CNFs can be obtained through low-intensity, hazard-free, processes from widely available biomass. Thus, this approach offers a unique shift in the range of opportunities to produce CNFs facilitating the integration of their use into the food supply chain, biomedical applications, and other regulatory-constrained applications.

NSD2 is a histone methyltransferase that modifies lysine 36 in histone H3 (H3K36), playing a central role in chromatin organization and transcriptional regulation. Oncogenic mutations, such as E1099K and T1150A in NSD2, have been associated with hyperactive methylation, but the molecular mechanisms underlying this gain of function remain poorly understood. In this study, we performed all-atom molecular dynamics simulations on models of NSD2 bound to the nucleosome for the wild type (WT), E1099K, T1150A, and the E1099K/T1150A double mutant. Analysis of MD simulations reveals that the global dynamics of the enzymes remain unaltered upon mutations. The time-lagged independent component analysis (tICA) and Markov state modeling uncovered fundamental differences in free-energy landscapes among the variants. The WT NSD2 exhibited energetically and kinetically unfavorable transitions between the macrostates along with extended enzyme-substrate distances. On the other hand, the mutant systems demonstrate reduced SAM-H3K36 distances with modified energy landscapes that facilitate transitions or favor prolonged occupancy of catalytically competent states. Importantly, the mutations reorganize the network of intramolecular contacts around the catalytic site, SAM-binding pocket, and histone-binding interface, optimizing substance engagement geometry. These findings demonstrate that oncogenic mutations achieve hyperactivity through strategic reorganization of conformational dynamics rather than simple destabilization, balancing local flexibility with global stability to enhance catalytic efficiency. Our results provide mechanistic insights into NSD2 dysregulation in cancer and establish a framework to develop allosteric inhibitors that target the enzymes dynamic landscape.

Building on a complex between a tubulin protofilament (PF) and a fragment of the Tau protein containing residues 169 to 367, we investigate the dynamics of the disordered elements of the system, namely the tubulin C-terminal tails (CTTs) and the Tau protein, using classical all-atom molecular dynamics simulations. Our results show that CTTs adopt a hook-like dynamic pattern on the bare PF while remaining highly mobile. The binding of Tau on the PF surface alters the dynamics of the I-CTTs in a sequence-dependent manner. While the repeat domains of Tau are mostly maintained on the PF by weak and strong binding patches with the tubulin cores, the Proline-Rich Region (PRR) relies on the wrapping phenomenon of I-CTTs to fuzzily stabilize its interaction with the PF. Our study thus provides a deep dive into the dynamic interplay between the Tau protein and the CTTs of microtubules, the latter being characterized extensively using a variety of disorder-adapted metrics. TOC Graphic O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=111 SRC="FIGDIR/small/721901v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@b3f985org.highwire.dtl.DTLVardef@1c2bf70org.highwire.dtl.DTLVardef@a66b95org.highwire.dtl.DTLVardef@1e138e0_HPS_FORMAT_FIGEXP M_FIG C_FIG

Magnetic treatment in tomato seeds (Solanum lycopersicum L.) has been studied as a biotechnological technique to induce a reduction in germination times and enhance plant development. However, the modified cellular mechanisms involved in the reduction of germination times or the improvement of development parameters are not yet clearly established. To explore a possible altered cellular mechanism, the effect of homogeneous static magnetic fields on the structure of the cyclic nucleotide-gated channel 6 (CNGC6), the modification in the organization of POPC lipids in the plasma membrane, and changes in calcium ion mobility were evaluated. For this purpose, coarse-grained molecular dynamics simulations were performed using the Martini 3 model in GROMACS, applying five different magnetic flux densities (0.000, 0.001, 0.010, 0.100, 1.000, and 10.000) T over 1 000 ns. The results showed an anisotropic effect in the longitudinal direction of the protein, which generated heterogeneous behavior among the chains of the homotetramer; this altered the conformation of the CNGC6 channel and modified the pore bottleneck. In contrast, no significant changes were observed in the conformational order of the POPC phospholipid chains. As a preliminary, single-replicate exploratory study, these results suggest that homogeneous static magnetic fields may induce specific structural modifications in the CNGC6 ion channel of Solanum lycopersicum L. without compromising the integrity of the lipid bilayer or the dynamics of ion transport within the analyzed timescale; these preliminary findings provide a molecular-level structural basis for future experimental and computational investigations of magnetic field effects on plant cyclic nucleotide-gated channels.

Surface functionalization of inorganic quantum dot nanoparticles is of great interest in the application of these materials toward a wide range of biological applications where membrane interactions are critical. The use of amphiphilic lipids to functionalize the surfaces of quantum dots represents a promising alternative to produce water-soluble and membrane-active materials with facile tuning of the quantum dots surface properties. Here, we demonstrate an experimental approach that yields lipid-coated quantum dots with highly tunable surface charge by controlling the concentration of cationic lipids during preparation. Through fluorescence-activated cell sorting assays, we show that these cationic lipid-coated quantum dots can enhance membrane interactions and increase membrane labeling density in live HEK293 cells. We further employed coarse-grained molecular dynamics simulations to model the lipid self-assembly process using an implicit solvent force field and subsequently model the adsorption of lipid-coated quantum dots to model membranes. Our simulations show that we can control the effective surface charge of lipid-coated quantum dots and influence the strength of adsorption to oppositely charged lipid membranes, a process that is mediated by the release of counterions at the quantum dot-membrane interface. This work supports the future development of biocompatible and water-soluble inorganic nanoparticles with highly tunable surfaces, and provides mechanistic insight into how different lipids can influence nanoparticle-membrane interactions at a molecular scale.

Intrinsic tryptophan fluorescence is widely used as a sensitive reporter of protein conformational dynamics, yet the molecular origin of its temperature-dependent modulation remains unclear. Here we investigate the conformational dynamics of Trp134 in bovine serum albumin (BSA) using molecular dynamics (MD) simulations, free-energy calculations based on umbrella sampling and WHAM, quantum mechanical (QM) calculations, and QM/MM approaches. MD simulations show that the global structure of BSA remains stable while temperature induces a gradual population shift from the Ia+ to the Ia- rotamer. The corresponding free-energy landscapes reveal that this shift arises from subtle changes in basin stability and transition barriers along the rotameric coordinate. In contrast, standalone QM calculations on isolated tryptophan predict different energetic trends, highlighting the sensitivity of rotamer stability to electronic-structure treatments and environmental effects. QM/MM calculations partially reconcile these differences by incorporating the protein environment. Together, these results suggest that temperature reshapes the rotamer free-energy landscape of Trp134, leading to population shifts that modulate intrinsic tryptophan fluorescence in proteins.

Orthoflavivirus, such as West Nile Virus (WNV), dengue virus (DENV) and ZIKA virus (ZIKV), are globally distributed pathogens that pose substantial threats to human health. Currently, there are still no effective antiviral drugs for WNV or ZIKV. Despite the availability of two licensed DENV vaccines, their use remains limited due to potential risks, highlighting an urgent need for antiviral drug development. The highly conserved orthoflavivirus protease NS2B/NS3 is required for viral replication, making it a promising anti-flavivirus target. A major challenge, however, is the highly charged active site of this enzyme, which requires charged chemical matters with low bioavailability. An alternative and more attractive strategy is to target potential allosteric sites or folding intermediate states of the protease. In this work, we employ the topology-based coarse-grained G[o] modeling to explore the coupled binding and folding pathways of WNV NS2B/NS3 protease and study the effects of the widely used experimental construct with a G4SG4 linker between NS2B and NS3 on stability and folding. Our results provide a holistic conformational landscape of the protease binding and folding, including several key intermediate states. We find that the presence of the G4SG4 linker alters the folding pathways and destabilizes the NS2B C-terminus. The latter is consistent with experimental observations that the G4SG4 linked protease has lower activity and adopts an open state without the substrate in crystal structures. Together, these findings provide for the first time a complete picture of the binding and folding of the NS2B/NS3 protease and identify important folding intermediate states that could be targeted for allosteric antiviral drug development. TOC Figure O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=157 SRC="FIGDIR/small/722635v1_ufig1.gif" ALT="Figure 1"> View larger version (40K): org.highwire.dtl.DTLVardef@163c356org.highwire.dtl.DTLVardef@ad7b35org.highwire.dtl.DTLVardef@173ed8aorg.highwire.dtl.DTLVardef@1f026bf_HPS_FORMAT_FIGEXP M_FIG C_FIG

We present a comparative study of the inner membrane of three Gram-positive bacterial strains, namely S. aureus, S. epidermidis and N. lacusekhoensis. A lipidomics study is used to obtain the lipid architecture and composition for S. epidermidis found in the skin microbiome and N. lacusekhoensis, an extremophile present in halophilic and alkophilic environments. Differences between the strains arise from both the lipid architecture and the cardiolipin content varying from 5% in S. aureus to 85% in N. lacusekhoensis. We develop coarse grained (CG) Martini-3 membrane models which reproduce structural properties such as membrane area, thickness, density distributions as well as ion-correlations with all-atom models. Inter-lipid correlations reveal a homogeneous distribution of lipids in the membranes despite the wide variation in lipid types and composition. Mechanical properties such as the area stretch modulus increased with cardiolipin content, however the bending modulus has a more complex dependence on membrane charge and lipid type. Using the CG models we evaluate the insertion free energies for four widely used antimicrobial molecules. Entry barriers for thymol and methylparaben arise from the charge density modulation at the membrane headgroups due to counterion condensation. The entry mechanisms of the antimicrobial peptide cecropin-melittin-15 (CM15) and the preservative molecule ethyl-lauroyl-arginate (ELAR) are found to be similar across all three strains. We also illustrate the manner in which the extremophilic strain, N. lacusekhoensis with its high cardiolipin content, modulates the partitioning kinetics of the antimicrobial molecule thymol with pH and salt. Our study reveals that membrane properties are largely conserved across the three model membranes. The molecular models and insights emerging from the present work should aid in the development of novel antimicrobials against Gram-positive strains.

"Mirror life", self-replicating organisms composed of non-natural-chirality biomacromolecules, presents a future threat with potentially global consequences. Consequently, there is strong agreement among experts that it should not be created. However, there is some disagreement over how effective existing medical countermeasures might prove against mirror bacteria in the event that they were created. Here, we leverage computational chemistry methods including docking and molecular dynamics to determine the likely binding efficacy of existing antibiotics against natural and mirror bacterial protein targets. We find that most existing antibiotics fail to bind to mirror bacterial protein targets, unlike their natural-chirality targets. This suggests altered binding of current medical countermeasures, which may impact the antimicrobial activity against mirror bacteria were the latter were created.